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PREFACE 


This  report  is  the  fourth  of  a set  of  five  presenting  the  results 
of  part  of  the  work  done  under  contract  number  F -30602-71-0-0312. 

The  report  is  delivered  to  RADC  in  accordance  with  item  A006  of 
the  Contract  Data  Requirement  List.  Sponsorship  and  technical 
direction  of  this  task  originated  in  the  Reliability  arri  Maintain- 
ability Engineering  Section  (A.  Coppola,  Chief) , Reliability  Branch 
(D.  Barber,  Chief) , within  the  Reliability  and  Compatibility  Division 
(J.  Naresky,  Chief)  of  the  Rome  Air  Development  Center.  Mr.  Anthony 
Coppola  was  the  Project  Engineer  who  was  technically  supported  by 
Mr.  Jerome  Klion. 


The  titles  of  the  reports  on  the  subject  "Reliability  Acceptance 
Sampling  Plans  Based  Upon  Prior  Distribution"  are  as  follows: 

Volume  I.  Introduction  and  Problem  Definition. 

Volune  II.  Risk  Criteria  and  Their  Interpretation. 

Volume  III.  Implications  and  Determination  of  the  Prior  Distribution. 

Volume  IV.  Design  of  Testing  Plans. 

Volume  V.  Sensitivity  Analyses. 


iii 


□ c 


ABSTRACT 


Numerical  and  graphical  procedures  for  the  design  of  reliability 
acceptance  sampling  plans  are  developed  in  this  report.  The  failure 
distribution  is  considered  to  be  exponential  with  parameter  6 and 
the  prior  distribution  of  0 is  taken  to  be  inverted  garrma.  These  pro- 
cedures are  applicable  when  accept/r eject  decisions  are  to  be  made 
regarding  a sequence  of  systems  or  lots.  For  specified  ccrrbinaticns 
of  producer-consumer  risks,  the  numerical  and  the  graphical  procedures 
can  be  used  for  the  design  of  a truncated  single  sample  plan  for  a 
system.  Other  single  sample  plans  can  be  obtained  by  using  the 
equivalence  relations  derived  in  this  report.  The  use  and  advantages 
of  the  graphical  procedure  are  illustrated  via  numerical  examples. 
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1.  INTRODUCTION 


This  report  presents  numerical  and  graphical  procedures  fcr  the 
design  of  reliability  acceptance  sampling  plans  for  the  case  when  the 
failure  distribution  is  exponential,  i.e. 

f (t| 6)  = 9 ^exp  (-t/0);  t > 0,  9 > 0 

where  t denotes  the  time  of  failure  and  9 is  the  mean  time  between  failures 
(MTBF) . The  parameter  9 is  considered  to  be  a random  variable  with  known 
prior  frequency  distribution  g (9 ) . This  prior  distribution  is  assured 
to  be  the  conjugate  inverted  garma  density 

g ( 0)  = yX9“(X+1)  exp (-y/0)/r (X) 7 y,X,0  > 0. 

The  successive  lots  or  systems  are  taken  to  have  different  (constant  but 
unknown)  01,  i = 1,2,3,...;  dram  independently  from  g(9). 

These  procedures  are  applicable  when  acaept/reject  decisions  cure  to 
be  made  regarding  a sequence  of  systems  or  a sequence  of  lots.  It  is 
assumed  that  the  decisions  are  based  on  data  subject  to  random  fluctuations, 
that  a definite  measure  of  loss  is  associated  with  each  of  the  two 
decisions  when  they  are  inappropriately  taken  and  that  each  lot  or  system 
is  sentenced  on  its  own  merit  without  regard  to  the  decisions  taken  about 
other  lots  or  systems. 

Following  four  types  of  single  sample  plans  are  considered: 

(i)  truncated  plan  with  replacement  (ii)  truncated  plan  without  replacement 
(iii)  censored  plan  with  replacement  and  (iv)  censored  plan  without  replacement. 
For  a lot,  which  consists  of  a large  number  of  components,  all  the  four  plans 
are  feasible.  For  a system,  which  can  be  repaired  as  new  upon  failure, 
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100%  Inspection  is  assured  and  plans  without  replacenent  are  infeasible. 
For  specified  combinations  of  produoer-oonsuner  risks,  the  numerical  and 
the  graphical  procedures  can  be  used  for  the  design  of  a truncated  single 
sanple  plan  for  a system.  Other  single  sanple  plans  can  be  obtained  by 
using  the  equivalence  relations  derived  in  this  report. 

Ihe  use  and  advantages  of  the  graphical  procedure  are  illustrated 
via  four  exanples.  Specifically,  the  examples  illustrate  the  design  of 
plans  and  the  exploration  of  the  design  region,  use  of  multiple  risk 
criteria,  design  with  engineering  constraints  and  design  with  partial 
information. 


2.  EQUIVALENCE  DELATIONS 


Plans  with  identical  P(a|0)  vs  0 curves  are  identical  i.e.  plans  with 
the  same  operating  characteristic  (O.C.)  are  equivalent.  Wte  now  determine 
the  conditions  under  which  the  different  types  of  single  sairple  plans 
have  the  sane  operating  characteristic  curve  as  the  truncated  plan  for 
a system.  That  is,  we  derive  the  equivalence  relationships  between  a 
single  sample  truncated  plan  and  the  other  types  of  plans. 


2.1  Plans  with  Replacement 


2.1.1.  Truncated  Single  Sanple  Plan  For  A 

Under  this  scheme  a single  system  is  placed  on  a life  test  of  duration  T. 
Without  loss  of  generality,  repair  is  assumed  to  be  instantaneous.  If  the 
observed  number  of  failures,  r,  is  less  than  or  equal  to  a value  r*  i.e. 
r ^ r*,  the  system  is  accepted.  Otherwise,  it  is  rejected.  Since  the 
failure  distribution  is  exponential , the  number  of  failures  in  fixed  time 
T has  a Poisson  distribution.  Then,  the  probability  of  acceptance  for  a 
given  0,  can  be  written  as: 


P(A|0)  = 


r* 

E 

r=0 


(1) 


■ 


2.1.2.  Truncated  Plan  For  A Lot 

In  this  schane  a sanple  of  n items  from  a lot  is  placed  on  a life  test 
of  duration  tQ.  Failed  items  are  immediately  replaced.  If  r ^ r' , the 
lot  is  accepted,  otherwise  it  is  rejected. 

The  failure  process  can  be  considered  as  a Poisson  process  with  rate 
(n/6) . Hence, 

r*  -nV0  r 

P(A|6,n)  =P(r^r,|0,n)  = E e (ntp/6)  (2) 

r=0  rl 

It  is  easy  to  see  that  the  O.C.  curves  given  by  Equations  (1)  and  (2) 
axe  identical  if  ntg  = T and  r 1 = r*.  Thus,  the  desired  plan  for  this 
situation  will  be  given  by  (tg=  T/n,  and  r'  = r*) . A suitable  value  of 
n is  chosen  to  make  the  sanple  representative  of  the  lot  and/or  from  cost 
considerations . 


2.1.3  Censored  Plan  For  A System 


A system  is  life  tested  until  r ' failures  occur.  If  the  total  test 

time  t'  > T',  the  systan  is  accepted;  otherwise  it  is  rejected.  If 

r' 

t. , t_ . . . t . represent  the  times  to  failure  then  t'  = E t.  las  a ganma 

i=l  1 

distribution  given  by. 


f(t'|0) 


T e (r'-l> 


er  r(r') 


w 


Hence, 


P(A|6) 


P(t'  > T'|0)  = / 
T' 


' r'-l 
e 6 (t')r  1 

er'  r(r') 


dt' 


If  t"  - t'/9,  then 


P(A|6) 


go  -f " r*«1 

* e z t"r  A 


/ 

T'/e 


Itr') 


dt"  = I 
r=0 


r'-l  -T'/B 


(T'/B)r 


(3) 


r! 


Equations  (3)  and  (1)  are  identical  if  r'  = r*  + 1 and  T'  = T.  Hence 
the  desired  plan  is  (T,r*  +1) . 


2.1.4  Censored  Plan  For  A Lot 

A sanple  of  n items  from  a lot  is  life  tested  until  r' <:  n failures 
occur.  If  the  total  test  time  t’  > T' , the  lot  is  accepted;  otherwise 


it  is  rejected. 

2t'  2 

It  is  shewn  in  Epstein  (1960)  that  — q—  X^i 
Hence, 


P (A  1 6)  = P(t*  > T'  1 6)  = P(—  > — ) 

e e 


= / 


where 


x = 


2T'  T(r')  2r 

e 

2t* 

0 


1 r'-l  -x/2 . 

T — x e dx. 
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e 

Ccnparing  Equations  (1)  and  (4) , the  required  plan  is  obtained  by  taking 


r'  ■ r*  + 1 and  T*  = T. 


W 


2.2.2  Censored  Plan  For  A Lot 

n items  from  a lot  are  life  tested  until  r'  £ n failures  occur.  Failed 
items  cure  not  replaced.  If  the  total  teeti  time  t'  > T' , the  lot  is 
accepted;  otherwise  it  is  rejected. 

2 

Epstein  (1960)  shows  that,  as  in  the  replacement  case,  2t'/6  ^ X2r«  • 
Hence  the  designed  plan  is  r'  = r*  + 1 and  T'  = T. 

Hie  relationships  between  the  above  plans  are  sunmarized  in  Table  1. 
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3.  RISK  CRITERIA 


Since  the  parameter  6 is  assured  to  have  a prior  density  g ( 0) , the 
producer's  and  the  consumer's  risks  may  be  quantified  in  various  ways. 

The  definitions  of  the  various  risks  are  new  given  in  terms  of  P (A  1 0) 
and  g(0).  For  the  interrelationships  and  interpretations  of  various  risks, 
the  reader  is  referred  to  Goel  and  Joglekar  (1976  b) 

3.1  Produoer's  Risks 
Classical : 

Hie  probability  of  rejecting  a system  with  0 equal  to  the  specified 
value  0q  is  equal  to  a,  i.e. 


a * P(R|60),  or  1 - a = P(A|0Q) 


(6) 


Average: 

The  probability  of  rejecting  a good  system  (0  £ 0g)  is  equal  to  a,  i.e. 
a = P (R 1 0 >,  0O) 


or 


1 - a-  / P(A|9)  9(6)de 

t r g(e)d0 
e0  0rt 


(7) 
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Posterior: 


Bib  probability  of  a rejected  system  being  good  is  equal  to  a*,  i.e. 


a*  « p ( 6 > 0 | R) 


or 


./{1-P(A|0)}  g (6)  d9 


»*  = ~ 


(8) 


1-  /P(A|  0)g(0)d0 

o 


Probability  of  Rejection: 

The  probability  that  a randomly  selected  system  will  be  rejected  is 
prespecified.  This  probability  represents  the  long  range  fraction  of  the 
submitted  systems  that  will  be  rejected  by  the  sanpling  plan  and  is  given 
by 


P (R)  = 1 - / P(A|0)  g(0)d6 
0 


(9) 


3.2  Consumer 1 s Risks 
Classical: 


The  probability  of  accepting  a system  with  minimum  acceptable  MTBF  0^ 
is  equal  to  6,  i.e. 


8 = P(A|01) 


(10) 
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Average: 

The  probability  of  accepting  a system  of  unacceptable  reliability 
(0  ^ 0^)  is  equal  to  B,  i.e. 


or 


B = p(a|0  < 0^ 


B = 


P(A|0)  g(0)d0 
jf1  g(0)d6 


(11) 


Posterior: 

The  probability  of  an  accepted  system  being  bad  (6  ^ 9^)  is  equal 
to  6*,  i.e. 


6*  = p(0  ^ e |a) 


or 


/XP ( A | 0)  g ( 0 ) d0 


7p(A|  0)  g (0)d0 
o 


(12) 


l 
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/ 

j 

Alternate  Posterior: 

To  provide  a better  control  on  the  distribution  of  0 in  the  accepted 
systems,  we  define  an  additional  posterior  risk,  8**,  which  specifies  the 
probability  that  an  accepted  system  will  have  an  MTBF  less  than  the 
specified  value  0Q.  Thus, 

8**  = P(9  < 0q|a) 
or 

0 

/ °P(a|0)  g(0)d0 

0**  = JL (13) 

/ 'P(A|0)  g(0)d© 

0 

The  various  risk  criteria  have  been  defined  above  in  terms  of  P(a|0) 
and  g (0 ) . Explicit  expressions  for  the  risks  for  all  the  four  types  of  single 
sample  plans  can  be  obtained  by  substituting  the  corresponding  expressions 
for  P (A|  0 ) derived  earlier.  Risk  expressions  for  the  truncated  plan  for  a 
system  sure  given  below. 
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NUMERICAL  PROCEDURE  FOR  THE  CESIQJ  OF  PLWS 


Plan  design  far  a specified  pair  of  producer-consumer  risks  is  achieved 
by  solving  the  corresponding  pair  of  equations.  Alternately,  one  may  only 
consider  the  design  of  truncated  plan  for  a system  and  obtain  the  remaining 
plans  by  using  the  equivalence  relations  given  in  Table  1. 

As  an  example,  if  the  design  criteria  are  P(R)  and  6*  then  given 
(y, \,  0^  0Q,  P (R) , 6*)  Equations  (17)  and  (20)  can  be  solved  to  obtain  T and 
r*.  Let  T*  = T/0Q,  Y*  = Y/9p,  K = 9,/^  and  <P  = 9/0 Q.  Then  Equations  (17) 
and  (20)  become 


P (R) = 1-  / 


e-T*/4>(TV(t,)r  }Y*X  (X+l)-Y*/<t>c 


1fK  { l*  e-'r*/<t>(Ty<l>)r}  Y*X  t-(X-H)-Y*/4>  d(j, 

Q r=0 r I r (X) 

f T e~TV^(TV»)r_  Y*X  r(*+l ^ 


f { I 
0 r=0 


Clearly,  given  O'*,*,  K,  P(R),0*)  the  values  of  T and  r*  can  be  numerically 
obtained,  if  the  solution  exists. 

The  numerical  procedure  consists  of  searching  for  those  values  of  T 
and  r*  in  the  design  region  that  satisfy  the  desired  risk  criteria.  Let 
us  call  the  consumer's  risks  6*,  8**,  S’  as  the  left  hand  risks  (LHR)  and 
the  producer's  risks  a*,  a , p(R)  as  the  right  hand  risks  (RHR) . It  is 
intuitively  clear  that  the  LHR  will  increase  with  increasing  r*  and  de- 


j 


f 


15 


creasing  T;  while  the  LHR  will  increase  with  decreasing  r*  and  increasing  T. 
Therefore,  the  region  of  feasible  plans  is  a V-shaped  area  bounded  by  lhr  and 
RHR  as  shown  in  Figure  1. 

In  order  to  search  for  the  plan  for  the  specified  prior  and  specified  risk 
criteria  etc.,  the  search  region  for  T and  r*  is  first  delineated.  We  then 
set  r*  equal  to  its  lower  limit  and  use  a binary  search  method  to  find  a 
Tr  such  that  the  plan  (Tr,  r*)  yields  a IiiR  closest  to  the  desired 
value,  the  closeness  being  specified  a-priori.  For  example,  we  can  require 

IHR  of  (T  ,r*)  - desired  LHR 

$.01 

Desired  LHR 

For  these  values  of  (T^ , r* ) , the  corresponding  RHR  is  canputed.  If  the 
canputed  value  is  less  than  or  equal  to  the  desired  value  then  the  plan  found 
is  the  desired  one.  Otherwise  we  increase  r*  by  1 unit  and  continue  the 
search  by  repeating  the  above  procedure.  The  search  is  continued  until  we 
find  a plan  or  we  do  not  find  one  in  the  delineated  region.  In  the  latter 
case  we  increase  the  size  of  the  region  and  proceed  as  discussed  above. 

A similar  procedure  is  followed  if  the  RHR  is  first  canputed  and  the 
test  for  meeting  the  LHR  criterion  is  satisfied  next.  The  plans  so 
obtained  will  not  be  exactly  the  same  in  general  but  will  approximately 
satisfy  the  desired  risks. 

A description  of  the  SUBROUTINE  DSP  for  designing  single  sarrple 
truncated  plans  using  the  above  procedure  is  given  in  Appendix  A. 
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Given  (Y*,X,  K,  a',  B'),  where  y*  « y/0q,  tlie  graphical  pro- 
cedure permits  the  computation  of  T*  and  r*  where  a'  and  6'  denote  the 
desired  producer's  and  the  consumer's  risk  respectively.  The  procedure 
consists  of  first  obtaining  contour  nomograms  by  drawing  contours  of 
constant  a'  and  8'  in  the  T*  - r*  plane.  The  contours  are  obtained 
as  follows.  In  the  T*  - r*  plane  a grid  is  constructed  with  r*  = 0(1)10 
and  T*  = 0(0.1)  9.  Far  given  X and  y*  and  K,  the  values  of  a'  and  B' 
are  evaluated  at  each  of  the  grid  points  using  SUBROUTINE  RISK,  given 
in  Appendix  B.  Finally,  contours  of  constant  a'  and  8'  are  drawn  by 
parabolic  interpolation.  These  contours  can  be  used  to  obtain  T*  and  r*, 
for  the  desired  combination  of  a'  and  8' . It  should  be  noted  that  a 
new  set  of  contours  will  be  required  for  each  pair  of  prior  parameters 
(X,y*) . However,  such  contours  can  be  easily  obtained  by  using  SUB- 
ROUTINE RISK  and  any  subroutine  for  drawing  contours  ttat  may  be  available 
on  a ccnputer  system. 

As  an  example,  of  the  design  of  a test  plan  using  the  graphical 
procedure,  consider  the  design  criteria  a = 0.05  and  8*  = 0.05. 

Let  the  prior  parameters  be  y = 1000  and  X = 0.4.  Let  0q  = 4000  and 
01  * 2000  such  that  the  discrimination  ratio  K = 2.  The  plots  of  constant 
a and  8*  are  shewn  in  Fig.  2.  For  a = 0.05  and  B*  = 0.05,  the  shaded 
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zone  in  Figure  2 represents  the  solution  region,  namely,  any  pair  of 
T*  and  r*  values  in  this  region  inplies  P(A|0  > eQ)  ^ 0.05  and 
p(6  $ el|A)  $ 0.05.  The  smallest  permissible  value  of  r*  is  1 and 
the  corresponding  value  of  T*  is  1.  Therefore,  the  designed  plan  is 
r*  - 1,  T « 4000.  Various  other  types  of  single  sample  plans  can  be 
ccrputad  from  Table  1 and  are  given  below. 

Truncated  Plan  for  a System:  r*  = 1,  T * 4000 

Censored  Plan  for  a System:  r*  = 2,  T * 4000 

Truncated  Plan  For  a Lot  With  Replacement:  n = 10,  r*  = 1,  tQ  = 400 

Censored  Plan  For  a Lot  With  Replacement:  r*  = 2.  T = 4000 

I Truncated  Plan  For  a Lot  Without  Replacement:  n = 10,  r*  = 1,  tg  = 421 

Censored  Plan  for  a Lot  Without  Replacement:  r*  * 2,  T = 4000 

The  use  and  advantages  of  this  procedure  are  now  illustrated  through  several 
numerical  examples. 
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An  advantage  of  the  graphical  procedure  is  that  it  permits  an  exploration 
of  the  T*  - r*  plane  near  the  solution  region.  Fig.  3 gives  the  contours 
of  constant  P(R)  and  3*  for  y*  = 0.125,  X » 0.2  and  K = 2.0.  If  the 
design  criteria  are  P(R)  $ 0.2  and  3*  •*  0.05,  then  the  shaded  region 
represents  the  solution  region  f ran  which  the  following  table  of  feasible  T* , 
r*  values  is  constructed. 


T* 

£>(rh 


1.06 


1.15 


1.42 

1.76 


1.65 

2.11 


2.10 


2.60 


The  notation  T*  ^ ^ represents  the  value  of  T*  which  satisfies  the 
specified  risk. 
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Corresponding  to  a given  value  of  r*,  the  feasible  values  of  T*  lie 
between  T* and  T*^R^.  When  T*  « ^(g*)'  ®*  = 0.05  P(R)  < 0.2 
while  T*  » T*  , inplies  that  P(R)  * 0.2  and  6*  < 0.05.  Ftor  other  feasible 

V*  / 

values  of  T*  in  the  feasible  region,  both  risks  satisfy  the  inequalities 
P(R)  N<  0.2  and  6*  * 0.05. 

Fran  the  viewpoint  of  testing  costs,  the  solution  with  mininun  T*  and 
r*  is  the  best.  For  the  present  exanple,  such  a solution  is  r*  = 2, 

T*  * 1.06.  Referring  to  Fig.  3,  it  is  of  interest  to  note  that  the  plan 
r*  = 1,  T*  = 0.7  very  nearly  satisfies  the  risks  with  over  33%  reduction 
in  test  time.  This  would  not  be  obvious  if  the  plans  were  designed  numerically 
or  presented  in  a tabular  form. 

5.2  Multiple  Risk  Criteria: 

The  graphical  procedure  is  easily  extended  to  the  specification  of 
nultiple  producer-consumer  risks.  In  Fig.  4,  the  example  of  Section  5.1 
is  reconsidered  with  the  added  consumer's  risk  B**.  Let  the  risk  criteria 
be  P(R)  4 0.2,  8*  ^ 0.05  and  8**  < 0.125.  These  criteria  ensure  that  at 
most  20%  of  the  systems  are  rejected,  at  most  5%  of  the  accepted  systems  have 
6 ^ 6^  and  at  most  12.5%  of  the  accepted  systems  have  6 ^ 6g. 

For  these  criteria  the  shaded  area  in  Figure  4 represents  the  solution 
region.  Ccnpared  to  the  solution  region  of  the  previous  example,  the  present 
solution  region  is  constrained  and  the  dominating  criteria  are  P(R)  and  8*** 
From  Figure  4,  the  designed  plan  for  this  case  is  r*  = 3,  T*  = 1.6. 
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Figure  4.  Graphical  Design  for  Multiple  Risk  criteria,  (P(R) 


The  graphical  procedure  also  indicates  the  feasibility  of  specified  risk 


criteria.  For  exanple,  if  0**  ^ 0.20,  then  this  criterion  can  not  be  met 
in  practice  for  the  desired  values  of  0*  and  P(R).  If  0**  > 0.20,  however,  the 
governing  risks  are  P(R)  and  3*  and  the  test  plan  derived  in  Section  5.1 
remains  unchanged. 

An  exanple  of  two  producer's  risks  is  given  in  Fig.  5.  The  criteria 
are  P(R)  0.2,  a $ 0.05  and  3*  ^ 0.05.  Solution  regions  for  (P(R),0*)  and  for 
(5,8*)  criteria  are  as  shown.  Clearly,  the  P(R),  8*  criterion  governs  and  the 
designed  plan  is  r*  - 2 and  T*  = 1.25. 


5.3  Design  With  Engineering  Constraints: 


At  times  engineering  considerations  will  require  that  testing  be  dis- 
continued after  a prefixed  time  or  that  only  a prespecified  n inter  of 
failures  be  allowed.  The  oontour  nanograms  can  be  used  to  design  an 
appropriate  plan  under  these  constraints.  The  procedure  is  illustrated  by 
reconsidering  the  exanple  of  Section  5.1,  and  referring  to  Figure  4. 


(a)  Let  us  assume  that  T*  i 1.  The  following  table  is  cons  true  ted  from  Fig.  4. 


T* 

P (R) 

8* 

in 

• 

o 

0.150 

0.080 

0.75 

0.205 

0.050 

1.0 

0.230 

0.030 

0.5 

0.105 

0.120 

0.75 

0.140 

0.085 

1.0 

0.180 

0.060 
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A reasonable  choice  of  T*,  8*  can  new  be  made.  For  exanple,  if  the 
desired  criteria  were  P(R)  £ 0.2,  8*  0.05  then  r*  = 1,  T*  = 0.75 

would  be  the  appropriate  plan. 


(b)  Let  us  now  assure  that  at  most  two  failures  are  permitted.  This  implies 
that  we  are  interested  in  a censored  plan  (2,T*) . The  corresponding 
truncated  plan  is  (1,T*) . The  following  table  is  constructed  from 
Fig.  4 to  determine  an  appropriate  T*. 


r* 

T* 

P(R) 

8* 

1 

0.5 

0.150 

0.080 

1 

0.75 

0.205 

0.050 

1 

1.0 

0.230 

0.030 

1 

1.25 

0.260 

0.018 

1 

1.5 

0.285 

0.011 

Again,  if  the  desired  risks  were  P(R)  ^ 0.20  and  8*  ^ 0.05,  the  plan 
is  r*  - 1,  T*  = 0.75,  the  corresponding  censored  plan  being  r*  = 2,  T*  = 0.75. 
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5.4  Design  With  Partial  Information: 

In  same  cases  only  one  of  the  tvro  risks  may  be  specified.  For  exaiqple, 
a producer  may  want  P(R)  ^ 0.2  but  the  consumer  is  not  certain  of  the 
consumer's  risk  and  vice-versa.  In  such  situations  a set  of  alternatives 
may  be  put  forth  as  follows. 


(a)  Suppose  P (R)  -$0.2  and  6*  is  unknown.  Then  the  following  table  can  be 
constructed  from  Figure  4 by  considering  the  points  of  intersection  of 
various  r*  values  and  the  P(R)  = 0.2  contour. 


r* 

T* 

8* 

8 

4.0 

0.026 

7 

3.53 

0.028 

6 

3.08 

0.030 

5 

2.60 

0.032 

4 

2.10 

0.035 

3 

1.65 

0.040 

2 

1.16 

0.044 

1 

0.70 

0.051 

(b)  Suppose  3*  .<  0.05  and  P (R)  is  unspecified.  The  following  table  can  be 
constructed  from  Figure  4 by  considering  the  points  of  intersection  of 
various  r*  values  and  the  3*  = 0.05  contour  . 


r* 

T* 

(R) 

6 

2.45 

0.166 

5 

2.10 

0.169 

4 

1.75 

0.173 

3 

1.40 

C.181 

2 

1.07 

0.190 

1 

0.72 

0.201 

Again,  a suitable  plan  can  be  selected  from  this  table. 


6.  CONCLUSIONS 


1.  Equivalence  relations  have  been  derived  in  this  report  between  the 
single  sanple  truncated,  censored,  with  replacement  and  without 
replacement  plans.  These  relations  permit  the  design  of  all  types 
of  single  sample  plans  cnce  the  truncated  plan  for  a system  is 
known. 

2.  A numerical  procedure  has  been  presented  for  the  design  of  single 
sample  truncated  plans,  vhen  the  parameter  0 of  the  exponential 
failure  distribution  has  an  inverted  ganma  prior.  This  procedure 
can  be  used  to  obtain  T and  r*  for  any  specified  producer's  and 
consumer's  risks  as  given  in  Equations  (14)  to  (21). 

3.  A graphical  procedure  for  the  design  of  truncated  single  sanple  plan 
for  a system  is  presented.  The  procedure  is  applicable  for  any 
specified  risk  criteria.  Other  single  sample  plans  can  be  obtained 
by  using  the  equivalence  relations. 

4.  Examples  are  presented  to  illustrate  the  advantages  of  the  graphical 
procedure.  Specifically,  the  examples  deal  with  the  design  and  ex- 
ploration of  the  solution  region,  specification  of  multiple  risk 
criteria,  design  under  engineering  constraints  and  design  with 
partial  information. 
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APPENDIX  A 


SUBROUTINE  DSP 


PURPOSE 

The  purpose  of  this  SUBROUTINE  is  to  search  the  testing  plan  (T,r*) 
in  a given  area  for  specified  parameters  of  the  inverted  gamma  prior 
distribution,  MXBF  9q,  discrimination  ratio  , and  the  desired  consumer's 
and  producer's  risks.  The  failure  time  distribution  is  assumed  to  be 
exponential. 

The  SUBROUTINE  also  gives  the  probabilities  P(0  4 0^  and  P(0  $ 0Q) 
USAGE 

CALL  DSP  0<L,KR,CX,CY,THO,DR,A/B,TIM,TRM/IR,IRU) 

DESCRIPTION  OF  PARAMETERS 


KL 

- the  consumer's  risk 

1 

for  6* 

2 

for  S' 

3 

for  e** 

KR 

- the  producer's  risk 

1 

for  a* 

2 

for  a 

3 

for  P (A) 

CX 

- the  specified  value  of 

the 

consumer's  risk 

CY 

- the  specified  value  of 

the 

producer's  risk 

THO 

- equals  0Q 
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DR  - discrimination  ratio,  Qq/^ 

A - scale  parameter  of  the  inverted  ganxna  distribution 

B - shape  parameter  of  the  inverted  ganna  distribution 

TIM  - the  lower  limit  for  searching  test  time 

TPM  - the  upper  limit  for  searching  test  time 

IR  - the  lower  limit  for  searching  acceptance  number 

IRLJ  - the  upper  limit  for  searching  acceptance  nuroer 

SUBROUTINES  REQUIRED 

INVG,  MDGAM,  UEKTST,  RISK 

LISTING 

A listing  of  the  SUBROUTINE  is  given  on L the  following  pages. 
METHOD 

The  search  method  for  obtaining  the  test  plan  (T,r*)  is  described 
in  Section  4 of  this  report. 


SUBROUTINE  DSP <KL t KR » CX  » CY  » THO » DR » A » B » TLH r TRM f IR r IRU ) 
DIMENSION  X(3)rY(3> 

EQUIVALENCE (X(1)fBS)f<X(2/f  QBQ )»<X<3)» B2S ) 

EQUIVALENCE ( Y (1)fAS)f(Y(2) fQAQ) f < Y(3) rPR) 

TRR=TRM 

IF(KR.EQ.3)  CY-l.-CY 

TH1-THO/DR 

ACC=»01 

CALL  INVG(AfBfTHOfPLTHO) 

CALL  INVG(AfBfTHIfPLTHI) 

TYPE  910fPLTH1fPLTH0 

910  FORMAT<///'  PR.  .LE.  TH1  = 'fF7.4f'  PR.  .LE.  THO=  'fF7.4> 
C***»  ACCEPT  WITHOUT  TESTING 
PA=1 . 

PR=1.-PA 

AS=0. 

BS=PLTH1 

QAQ-O. 

QBQ=1 » 

B2S=PLTH0 

IF  (X(KL).GE.CX)  GO  TO  90 
IF  (YtKR) .GT.CY)  GO  TO  90 
TYPE  930 

930  FORMAT ( ' ******  ACCEPT  WITHOUT  TESTING  ') 

TYPE  951fASfBSfOAQfQBQfPAfB2S 
RETURN 

C90****  REJECT  WITHOUT  TESTING 
90  PA=0. 

PR=1 • -PA 
AS=1 . -PLTHO 
BS=0. 

QAQ=! . 

QBQ=0. 

B2S=0. 

IF  (X(KL) .GT.CX)  GO  TO  100 

IF  (Y<KR) .GT.CY)  GO  TO  100 

TYPE  920 

970  format ( * ****  REJECT  WITHOUT  TESTING  ) 

TYPE  951 fASfBSfQAQfQBQfPAfB2S 
RETURN 

C *****STAR  SEARCHING  SOLUTION 
100  TRM=TRR 
T=TLM 

CALL  RISK(T»IRfTH0»DRfA»B»AS»BSfQAQfQBQ»PA»B2S) 

PR= 1 . -FA 

IF  (CY.GE.Y(KR))  GO  TO  12 

IF  (CX.GE.X(KL) ) GO  TO  401 

IF  (IR.GE.IRU)  GO  TO  250 
IR=IR+1 

cn  to  ioo 

IF  ( CX . GE » X ( KL ) ) GO  TO  101 
T=TRM 

CALL  RISK  < T f IRf  THOfDR  f Ar  Bf  AS . BSfQAQfGBQfPAfB2S ) 

PR=1 .-PA 

IF  ( X (KL ) . GE • CX ) GO  TO  200 
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IF  <CY.GE.Y<KR>)  GO  TO  10 
GO  TO  40 

101  TRM=TLM 
TLM=0. 

GO  TO  10 

401  TRM=TLM 
TLM=0. 

GO  TO  40 

40  T=(TLM+TRM)/2. 

CALL  RISK  (Tf IR  f THOf  DRf A fBf  AS f BS  fQAGf  QBQ  fPA r B2S) 

PR=1.-PA 

IF  <CY »GE»Y<KR) ) GO  TO  20 
IF  (CX.GE.X(KL) > GO  TO  44 
TLM=T 

IF  (IR.GE.IRU)  GO  TO  250 
IR=IR+1 
GO  TO  100 
44  TRM=T 
GO  TO  40 

20  IF(CX.GE.X(KL) ) GO  TO  102 
TLM=T 
GO  TO  40 

102  TRM=T 

10  T=(TLM+TRM)/2. 

CALL  RISK  ( T f IRf  THO f DR  f Af  B f AS  f BS  f QAQ t QBQf  PA  fB2S ) 

PR=1 • -PA 

IF  (CX.GE.X(KL)  ) GO  TO  111 

TLM=T 

GO  TO  10 

111  IF  ( (ABS(CX-X<KL) ) )/CX  .LT.  ACC)  GO  TO  800 
TRM=T 
GO  TO  10 
250  TYPE  255 

255  FORMAT ( ' ****  NO  SOLUTION  IN  THIS  AREA » ' / 

*'  INCREASE  NUMBER  OF  ACCEPTANCE'/ 

* /'  THE  STOPPING  POINT  IS') 

TYPE  900 » T ? IR 

TYPE  951 »ASfBS»QAQfQBQfPA»B2S 
RETURN 

200  TYPE  210 

210  FORMAT < * ****  NO  SOLUTION  IN  THIS  AREA'/ 

*/'  INCREASE  TESTING  TIME'/ 

*/'  ****  THE ,J5J0PPING  POINT  IS') 

TYPE  900»TfIR^' 

TYPE  951 »AS»BS»QAQf GBQ»PAfB2S 
RETURN 

800  TYPE  810 

810  FORMAT  ('  *****  TESTING  PLANE  IS  THE  FOLLOWING') 

TYPE  900 fT » IR 

TYPE  951 f ASfBSfQAQfQBQfPAfB2S 
RETURN 

951  FORMAT < ' A*= 'fFA.4f2Xf/B*='fF6.4f2Xf ' A-= ' f F6 . 4 f 2X f 'B-= ' f F6 . 4 f 2X » 
A ' PA= ' f F6 . 4 f 3X  f ' B2S= ' f FA . 4 ) 

900  FORMAT ( ' T= ' fF12 . 3 f 3X f ' R*= ' f 13 f 3X f ' X= ' fF8 . 5 f 3X f ' Y= ' fF8 . 5 ) 

RETURN 

END 
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Given  a = 3,  y = 300,  eQ  = 100,  DR  = 2,  and  let  the  risk  criteria  be 
P(A)  ^ 0.8  and  8*  ^ 0.025.  We  search  for  the  testing  plan  for  test 
time  between  0.005  and  150;  and  for  acceptance  nuttoer  between  0 and  3. 

Celling  Program  and  Output 

CALL  DSP ( 1 ,3t »02S» . 8 » 100 . f 2 . >300. * 3 ♦ f . 005 » ISO . f0f5> 

STOP 

END 


MAIN. 

fortran:  dsp 
DSP 

LOADING 


2 DSP  SK  SORE 
EXECUTION 


PR.  . LE . Till*  0 . 0620  PR.  .LE.  THU*  0.1232 
*.!:**  NO  SOLUTION  IN  THIS  AREA 

INCREASE  TESTING  TIME 

**#*  THE  STOPPING  POINT  IS 

r 150.  COO  R*:=  3 x = 

ax  .121.3  BS-*.  .0351*  A ■••••■  .0209  .5114  PA*  .8999  B2S* 


Thus  we  have  P(  0 $ 6^  = 0.0620  and  P(  ® < 0Q)  * 0.4232.  However, 
the  output  shows  no  feasible  testing  plan  in  the  specified  region  for  T and  r 
and  suggests  increasing  the  testing  time.  So  we  increase  the  vpper  limit 
of  test  time  to  300.  The  output  now  is  as  follows 

CALL  DSP ( 1 r 3 » . 025 » .8> 100. r2.  >300. >3. >150. *300. 1 0*3> 

STOP 

END 


exec  dsp  • rel  t ccx . f 4 » a ♦ rel » ririew  ♦ f el 

fortran:  ccx 

MAIN. 

LOADING 

CCX  3K  CORE 
EXECUTION 


PR.  .LE.  TH1-  0.0620  PR.  .LE.  THQ~  0.4232 
*****  TESTING  PLANE  IS  THE  FOLLOWING 
T=  190.652  R*-  3 X= 

A*=  .1499  E<*=  .0248  A~~  .0426  B-~  .3343  PA  .8362 

STOP 

END  OF  EXECUTION 

CPU  TIME:  0.65  ELAPSED  TIME:  4.63 
EXIT 


Thus  we  get  the  desired  plan  to  be 

T = 190.654,  and  r*  = 3. 


37 


I 


f 


APPENDIX  B 
SUBROUTINE  RISK 

PURPOSE 

This  SUBROUTINE  oonputts  the  risks  a*,  6*/  g**,  a,  8,  P(A)  for  given 
test  time  T,  acceptance  muter  r*,  specified  MTBF  0Q/  discrimination  ratio  K 
and  the  given  parameters  of  the  inverted  gaima  distribution.  The  failure 
time  distribution  is  assumed  to  be  exponential. 


USAGE 


CALL 


RISK  (T,IR,TH0,DR,A,B,AS,BS,QAB,QBQ,PA,R2S) 


DESCRIPTION  OF  PARAMETERS 


T 

IR 

TH0 

DR 

A 

B 

AS 

BS 

QAQ 

QBQ 

PA 

B2S 


- test  time 

- acceptance  number 

" 90 

- discrimination  ratio,  equals,  9^/9^ 

- scale  parameter  of  the  inverted  garrma  distribution 

- shape  parameter  of  the  inverted  gamma  distribution 
~ posterior  producer's  risk,  a* 

- posterior  consumer's  risk,  B* 

- average  producer's  risk,  a 

- average  consumer's  risk,  B 

- unconditioned  probability  of  acceptance,  P(A)  = 1 - P (R) 

- alternate  posterior  consumer's  risk,  8** 
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SUBROOTINES  REQUIRED 

<?tf«A,ZNVG,  MXJAM,  UERTST 

\ 

LISTING 

A listing  of  the  SUBROUTINE  is  given  on  the  following  page. 
FCTPD 

The  various  risks  are  computed  based  on  the  expressions  given  in 
Section  3 of  this  report. 
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SUBROUTINE  RISK(Tr  IRrTHOf  DR»A>Br  ASrBS»QAQ»QBQrF‘A»B2S> 
DIMENSION  PI (99 ) 

IK=IR+1 

TH1=THG/DR 

CALL  INVG ( A » B r THO  »PLTHO) 

CALL  INVG(A»BfTH1»PLTH1 ' 

AB=A**B 

AT=A+T 

ATDTHO= A T /THO 

ATDTH1=AT/TH1 

ADAT=A/AT 

TDAT=T/AT 

ADATB=ADAT**B 

PI<1>=1. 

DO  10  1=2 » IK 
AI=FL0AT(I-1) 

10  PKI>=PI<I  -1)*TDAT*<B+AI~1.  >/AI 
PAA=0 
PB=0 . 

PB0--0 • 

DO  20  1IX=1»IK 
JJX=IIX~1 
TJX=FLOAT( JJX) 

W3=T  JX+B 

CALL  MDGAM ( ATDTH1 »U3»GITH1»IER) 

CALL  MDGAM < ATDTHO , U3 » GI INOfIER) 

PAX=PI ( I IX  > 

PJ1=PAX*<1.-GITH1 ) 

PJO=PAX*( 1 .-GITHO) 

PB-PB+PJ1 
PBG=PBO+PJO 
PAA=PAAiPI < 1 IX) 

20  CONTINUE 
PB=PB*ADATB 
PBO -PBO#ADA I B 
PA=PAA*ADATB 
BS=PB/PA 
PR-1 .-PA 

AS=1  »-<PLTHO--PBO)/PR 
QAQ= 1 . - < PA-PBO ) / ( 1 . -PLTHO ) 

QBQ=PB/PLTH1 

B2S=PB0/PA 

RETURN 

END 


Example: 

Given  X = 3,  y = 300,  0 Q = 100,  DR  = 2 and  the  testing  plan,  T = 190.654, 
r*  » 3,  we  want  to  obtain  the  various  producer ' s and  the  consumer's  risks. 


Calling  Program  and  Output 


CALL  R I CK  < 1 VO . 654 > 3 > 1 0 0 . > 2.  >300.  >3.  > AS  > BS  > QAG  r QBQ  > PA  > B2S ) 
TYPE  1 0 > AS  > BS  > Cl  AW  > QBW  > PA  > B2S 
10  rORMAT ( ' ' >6F10 . 4 ; 

STOP 

END 


CCX  2K  CORE 
EXECUTION 

0.1 1?y  0.0240  0.0426  0.3343  0.8362  0.3396 

STOP 

EMU  OP  EXECUTION 

CPU  IIMEJ  0.08  ELAPSED  (IMEJ  3.73 
EXIT 


Thus,  for  the  specified  prior  distribution,  6 q,  DR  and  the  test  plan, 
the  various  risks  as  ocrputed  above  are.* 

o*  = 0.1499,  0*  = 0.0248,  a = 0.0426,  B = 0.3343,  P(A)  = 0.8362, 

8**  m 0.3396. 
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APPENDIX  C 


DESCRIPTION  OF  PARAMETERS 

XX  - the  argument  for  the  garma  function 
GX  - the  resultant  gamma  function  value 
IER  - resultant  error  code  where 


IER  = 0 no  error 

IER  =1  XX  is  within  0.000001  of  being  a negative  integer 

IER  -2  XX  is  greater  than  57,  overflow,  GX  set  to  1.0E75 


mmmmmgmmmm 


SUBROUTINE  GMMMA (XXrGXxIER) 
IF(XX-57. >6,6,4 
4 IER=2 
GX=1 . E75 
RETURN 
6 X=XX 

ERR=1 ♦ 0E~6 
IER=0 
GX=1 .0 

IF ( X-2 .0)50»50fi5 
10  I F ( X-2 « 0 ) 1 1 0 y 1 1 0 y 1 5 
15  X=X~1.0 
GX=GX*X 
GO  TO  10 

50  IF ( X-l ♦ 0 ) 60 » 120 p 1 10 


SEE  IF  X IS  NEAR  NEGATIVE  INTEGER  OR  ZERO 

60  IF ( X-ERR ) 62  y 62  y 80 
62  Y-FL.OAT ( INT(X ) )-X 

I F ( A B S < Y ) - E R R) 130, 130, 6 4 
64  IF  < 1 * 0-Y--ERR ) 1 30 , 130  r 70 

X NOT  NEAR  A NEGATIVE  INTEGER  OR  ZERO 


70  I F ( X - 1 ♦ 0 ) 80  y SO  y 1 1 0 
80  GX-GX/X 
X-XH  « 0 
GO  TO  70 
110  Y=X  -1 ♦ 0 

G Y = 1 . 0 f Y * ( - 0 * 5 7 7 1 0 1 7 +Y*< + 0*985854 0 + Y * < - 0 * 8 7 6 4 2 1 8 + Y Y ( 
1 Y* ( -0 . 5684729+Y* ( +0 . 2548205+Y* < -0  * 051 49930 ))>)))) 
GX=GX*GY 
120  RETURN 
130  IER'1 
RETURN 
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APPENDIX  D 
SUBROUTINE  INVG 

PURPOSE 

This  SUBROUTINE  ccnputes  the  ciinulative  distribution  function  of  the 
inverted  gamma  density  f(0;A,B)  for  a given  argument  0 v.ith  scale  parameter 
A,  and  slepe  parameter  B. 

USAGE 

CALL  INVG  (A,B,TH,P) 

DESCRIPTION  OF  PARAMETERS 

A - scale  parameter  of  the  inverted  gamma  distribution 

B - shape  parameter  of  the  inverted  gamma  distribution 

TH  - 0 

P - P(0  ^ 0) 


SUBROUTINES  REQUIRED 
MX3AM 


LISTING 

This  SUBROUTINE  is  taken  from  the  International  Mathematical  and 
Statistical  Libraries. 


SUBROLJT I NE  INVG ( A y B y Tile  P ) 
X=A/TH 

CALL  M D G A M ( X y B v P » I E R ) 

P=l.-P 

RETURN 

END 
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APPENDIX  E 
SUBROUTINE  MDGAM 


PURPOSE 

MDGAM  evaluates  the  probability  that  a randan  variable  from  a gamma 
distribution  having  the  parameter  P is  less  than  or  equal  to  X.  Unis  MDGAM 
ocrrputes  the  incomplete  gamma  ratio 

icx,p)  - /X  e ttp"1  dt/r(p) 

0 

where  r CP)  is  the  ccnplete  gamma  function. 

usact; 

CALL  MDGAM  (X,P,PROB,IER) 


PARAMETERS 
X 
P 

PRDB 

IER 


value  to  which  gamra  is  to  be  integrated 
input  gamma  parameter 

output  probability  = integral  of  gamma  (P)  to  X 

error  indicator 

Terminal  error  = 128  + N 

N = 1 indicates  X is  less  than  zero 

N - 2 indicates  P is  less  than  or  equal  to  zero 


; 

LISTING 

The  listing  is  given  on  the  following  page. 

This  SUBROUTINE  is  taken  from  the  International  Mathematical  and 
Statistical  Libraries. 
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SUBROUTINE  MDGAM  <X»P»PROBf IER) 


5 


C 


c 


c 

c 


DIMENSION  V(6)»V1(6> 

EQUIVALENCE  (V(3)»V1(1)) 

TEST  X AND  P 

PROB  = 0.0 

IF  <X  .GE.  0.0)  GO  TO  5 
IER  = 129 
GO  TO  9000 

5 IF  <P  .GT.  0.0)  GO  TO  10 

IER  = 130 

GO  TO  9000 

10  IER  = 0 

IF  <X  .EQ.  0.0)  GO  TO  9005 

DEFINE  LOG-GAMMA  AND  INITIALIZE 

PNLG  = ALGAMA(P) 

CNT  = P * ALOG(X) 

YCNT  = X + PNLG 

IF  < ( CNT-YCNT ) .GT.  -174.673)  GO  TO  15 
AX  = 0.0 
GO  TO  20 

15  AX  = EXP (CNT-YCNT) 

20  BIG  = 1.E35 
CUT  = l.E-8 

CHOOSE  ALGORITHMIC  METHOD 
IF  <<X  .LE.  1.0)  .OR.  (X  .LT.  P ))  GO  TO  40 

CONTINUED  FRACTION  EXPANSION 

Y = 1.0  - P 

Z - X + Y + 1.0 
CNT  =0.0 
V(l>  = 1.0 

V<2)  = X 

V( 3 ) = X + 1.0 

V( 4 ) = Z * X 

PROB  = V(3)/V(4) 

25  CNT  = CNT  + 1.0 

Y = Y + 1.0 
Z = Z + 2.0 
YCNT  = Y * CNT 

V( 5 ) = Vl(l)  * Z - V(l)  * YCNT 
V<6)  = VI (2)  * Z - V(2 ) * YCNT 
IF  ( V ( 6 ) .EQ.  0.0)  GO  TO  50 
RATIO  = V(5)/V(6) 

REDUC  = ABS  < PROB-RATIO ) 

IF  (REDUC  .GT.  CUT)  GO  TO  30 
IF  (REDUC  .LE.  RATIO*CUT)  GO  TO  35 
30  PROB  = RATIO 
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GO  TO  50 

35  PROP  = 1.0  - PROP  * AX 
GO  TO  9005 

: SERIES  EXPANSION 

40  RATIO  = P 
CNT  =1.0 
PROP  =1.0 

45  RATIO  = RATIO  + 1.0 
CNT  = CNT  * X/RATIO 
PROP  = PROP  + CNT 
IF  (CNT  .GT.  CUT)  GO  TO  45 
PROP  = PROP  * AX/P 
GO  TO  9005 
50  DO  55  1 = 1,4 

0<I)  = OKI) 

55  CONTINUE 

IF  (APS (0(5) ) .LT.  PIG)  GO  TO  25 

> SCALE  TERMS  DOWN  TO  PREVENT  OVERFLOW 

DO  60  1=1.4 

V(I)  = V( I )/PIG 
60  CONTINUE 
GO  TO  25 
9000  CONTINUE 

CALL  UERTST  ( IER . 6HMDGAM  ) 

9005  RETURN 
END 

FUNCTION  ALGAMA(P) 

CALL  GMMMA(P.GP.IER) 

ALGAMA=ALOG( GP ) 

RETURN 

END 


APPENDIX  F 


SUBROUTINE  UERTST 


PURPOSE 

To  print  a message  reflecting  the  error  detected  by  a SUBROUTINE 
USAGE 

CALL  UERTST  (IER,NAME) 

PARAMETERS 

IER  - error  parameter.  Type  + N where 

type  = 128  implies  terminal  error 
64  implies  warning  with  fix 
32  implies  warning 

N = error  code  relevant  to  calling  routine 
NAME  - input  vector  containing  the  name  of  the  calling  routine 
as  a six  character  literal  string. 

LISTING 

Given  an  the  following  page. 


1 


19 


S U B R 0 U TIN E:  IJ E R T 3T  <1 £ R • I A .11' 

) 

BIMENGITN  11  TT 

A , 

TBIT ( : ' 

INTEGER  NAME  ( 6 '/ 

INTEGER  UARDU,, 

:-r  v 

: • : i 1 • r D v 'i  R 

EQUIVAL.ENCE  r TT-.'i 

DATA  ( ( ITYT  \ 1 ' 

:<  w 

v ■< 

ANN ) • T i T '»  flJ.iR 

/ 

. WAI.'.' 

t ' 

■ warn  * 

' J 

IT  ?:r 

/ “T  r ^ < • 

; • ' • ' > 

'NGN  ' 

r»rrr  t tvi 

• . » ..  i • * . • . * • 

DATA  < I D L T ( I ) > I i r 1 ) / 32 * 6 
DATA  PRINT R / 6 ' 

\ V .1. 

n A • 

v un.'/ 

I£R2=IER 

IF  ( 1 ! : R 2 C-r  WARN  ' GO  TP 

r- 

in.;  .:rr  T": 

IER1 - 4 
GO  TO  20 

IF  ( IER:.?  < !..T . TERM' 

IER1  3 
GO  TO  20 

IF  (IER2  .LT.  WARD 

1 0 

w.m  Ml  M • ; . 

IER 1-2 

GO  10  20 

j.  : 

1ER1  1 

L. N ru » iC  * ' 

IER2  IER2  I BIT  1 1.  •: ! > 

rnTNi  e:  re  : 

WRITE  ( PR I NT R ' > > I TiT  < I • 
FORMAT  ( ' **  * INS!  'Dl 
K ' (IER  '.13. 

in 

, i ; , 

RETURN 

